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Quantitative dose-response modeling is 
an important contributor to cancer risk 
assessment; hence, it is a major factor 
in cancer risk management and the reg¬ 
ulatory process. The dose-response 
models now used in the regulatory 
process are overly simplistic, probabi¬ 
listic representations of highly complex 
biological phenomena; tiiiese models 
are not biological models. Several of 
these simplistic models provide similar 
fits to the high-dose data generated in 
chronic anii^ bioassays but provide 
dissimilar projections of risk at the 
lower doses of interest to man. Figure 1 
shows model extrapolations. 

TIk possibilities for the low-dose be¬ 
havior of a simplistic model can be so 
independent of the fit of that model to 
the experimental dau that an upper 
confidence limit, or upper bound, on 
the risk at a low dose can be orders of 
magnitude larger than the fitted model 
value. Figure 2 shows that the upper 
bounds are not responsive to the experi¬ 
mental data, and Figure 3 shows that 
experimental outcomes with vastly dif¬ 
ferent observed dose-tesponse patterns 
all have essentially the same upper 
bounds. 

The potency measures, such as unit 
risks and relative risks, cited by the 
regulatory agencies are based on upper 
bounds a^ not on fitted model values. 
These measurements do not differenti¬ 
ate between carcinogens on the basis of 
available experimental data about the 
Stapes of the dose-response relation¬ 
ship. The inability to differentiate be¬ 
tween risks is a serious shortcoming 
when all alternatives involve some risk. 
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In Older to obtain more useftil quanti¬ 
tative dose-response assessments, the 
plethora of simplistic models must be 
replaced by more biolc^cally reflective 
dose-response models that utilize the 
available scientific information. New 
dose-response modeling techniques 
can incorporate representations of the 
exposure in terms of dose scales based 
on cell turnover rates, repair processes, 
immune system responses, physio¬ 
logical and pharmacokinetic models of 
the absorption, delivery, metabolism, 
and elimination of chemicals. 

These new techniques can incorpo¬ 
rate distributions of individual back¬ 
ground exp^ure levels and individual 
susceptibilities to low levels of a chemi¬ 
cal. litey also can incorporate such fac¬ 
tors as low-dose lirtearity on a biologi¬ 
cally relevant scale; age-dependent 
changes in the number of target cells 
and time-dependent effects of cell pro¬ 


liferation on the number of intermedi¬ 
ate cells in a multistage process; and 
dose levels, susceptibilities, and back¬ 
ground exposures that are not necessar¬ 
ily constant over time. These new, 
more biologically reflective models are 
now described as state-of-the-art, and 
additional research should lead to even 
better dose-response models. 


The dose level should be expressed 
on a biologically relevant scale. For ex¬ 
ample, Figure 4 shows that the dose 
level conesponding to an inhalation ex¬ 
posure can be expressed in several 
terms; the concentration of the chemi¬ 
cal in the air inhaled; the total amount 
of the chemical inhaled over time; the 
amount of the chemical or its active 
metabolite reaching the target cells of a 
particular tissue; and the amount of 
chemical or metabolite interaction with 


DNA or the amount of cell damage that 
escapes repair. 

The dose scale can be expressed as 
the administered dose, the intermediate 
dose, the dose delivered to the target 
site, or the biologically effective dose. 
The biologically effective dose reflects 
not only the amount of the chemical or 
metabolite delivered to the target site 
but also the net generation of cancer- 
related activity, including cell turnover, 
DNA repair, and immune system re¬ 
sponses. 

There are advantages to studying bio¬ 
logical response. Using a biologically 
reflective dose instead of the adminis¬ 
tered dose, the modeler can incorporate 
known biology, physiology, and phar¬ 
macokinetics. This reduces the uncer¬ 
tainty and the number of unknowns in 
the dose-response model and aids ex¬ 
trapolations from high doses to low 
doses, from species to species, from 
one exposure route to another, and from 
one exposure time frame to another. 

If there is insufficient scientific 
knowledge to establish the complete 
transformation from administered dose 
to a biologically effective dose, then us¬ 
ing even a partial transformation or 
best approximation based on what is 
known should lead to a better risk as¬ 
sessment than using only the adminis¬ 
tered dose. 

Individual susceptibilities 

Individuals vary in their responses to 
a particular exposure to a potentially 
ha^dous chemical. Individual varia¬ 
tion is influenced by many factors, in¬ 
cluding variation in background expo¬ 
sures, genetic traits, preexisting 
diseases, and behavioral traits (7). An 
extension of the linear multistage model 
utilizes the biologically effective dose 
and reflects individual variations in 
background exposure and biological 
susceptibility to low levels of a chemi¬ 
cal. An analogous extension of a more 
biologically motivated model that em¬ 
phasizes the multistage carcinogenic 
process, age-dependence, and cell pro¬ 
liferation is described below. 

The background exposure for a pop¬ 
ulation of individuals, Dg, is a random 
variable; that is, the level of exposure 
from all sources, other than the one 
source being regulated, is potentially 
different for different members of the 
population. Each individual has his 
own particular value dg for the back¬ 
ground dose Do- 

An animal'^ background exposure is 
dependent upon such factors as envi¬ 
ronment and diet. A human's do can 
depend on these same factors as well as 
occupation, lifestyle, and hygiene. 

Similarly, because there is often in¬ 
terindividual variation in biological fac¬ 
tors such as metabolic rates and enzyme 
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levels, the biological susceptibility S to 
a particular chemical varies in a popu¬ 
lation Irom one individual to another. 
For a population, 5 is a random varia¬ 
ble, ai^ each individual has a particular 
value s for the susceptibility 5. 

An Individualized Response Model 
utilizes the biologically effective dose 
as the dose scale for dose-response 
modeling and reflects the interindivi¬ 
dual variation in the susceptibility S and 
the background dose Do- A description 
of how the Individualized Response 
Model works is provided in the box. 
Figure S provides an overview of the 
components and structure of the model. 

A reasonable approximation for the 
distribution of the background dose Do 
in the population may often be the 
lognormal distribution. There has been 
considerable historic success in using 
the lognormal distribution to represent 
the distributions of a wide variety of 
exposures and biological phenomena. 
Lognormal distributions agree with 
Do > 0. In addition, lognor^ distri¬ 
butions can be very skewed to match 
chemical carcinogens for which Do is 
usually veiy near zero. On the other 
hand, lognormal distributions can be 
nearly symmetrical and more like the 
normal distribution if Do is frequently 
moderately large. Similarly, the distri¬ 
bution of the susceptibility 5 may be 
approximated by a lognormal distribu¬ 
tion. 

The biological phenomena repre¬ 
sented by susceptibility and the way in 
which susceptibility enters into the car¬ 
cinogenic process will generally be dif- 
fetent for different chemicals. Some¬ 
times susceptibility can be expressed 
through a susceptibility frontier. The 
delivered dose corresponding to the 
transition from substantially lowered 
carcinogenic effectiveness to the dose 
region where the predominant mecha¬ 
nism resisting or suppressing carcino¬ 
genesis is overwhelm^ is the suscepti¬ 
bility frontier. 

For a populatioa of individuals, die 
susceptibili^ frontier F is a nmdom 
variable. E^h individual has his or her 
own particular vahie/for the suscepd- 
bility frontier E Genetic differences in 
meuibolic activation and detoxification, 
tissue sensitivity, repair efficiency, and 
immune system responses contribute to 
the individual variation in the value of 
the susceptibility frontier (2, 3). 

A susceptibiliQr frontier is different 
from a tolerance. If the delivered dose 
becomes greater than an individual^ 
toleiance, then cancer is a certainty; 
however, as the delivered dose becomes 
greater titan the individuals susceptibil¬ 
ity frontier, then the individuals cancer 
probability increases more rapidly but 
doesn’t immediately become a cer¬ 
tainty. 
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The concept of a susceptibility fron¬ 
tier is also different from that of a 
threshold. A threshold denotes a dose 
level below which absolutely no carcin¬ 
ogenesis is possible, whereas a suscep¬ 
tibility frontier indicates only a transi¬ 
tion l^m the normal functioning of the 
bodyS carcinogenic defense mecha¬ 
nisms to dose levels overwhelming 
those mechanisms. Also, a threshold is 
the same for every member of the pop¬ 
ulation, whereas ok susceptibility fron¬ 
tier differs for each individual. 

The concept of a susceptibility fron¬ 
tier can be illustrated in simple terms as 
Mows: Suppose that for a particular 
individual a comical exposure causes a 
particular number of specific moleculre 
to reach and invade a target cell. Also, 
supper that the target cells have, for 
simplicity, a single cancer defense 
mechanism relative to these invaders— 
say a detoxification process—and that 
the cell begins with 20 opportunities to 
carry out the detoxification of an invad¬ 
ing molecule. However, at each oppor¬ 
tunity the probability of detoxifying an 
inva^r is, for example, only 0.9 (not 
1.0), and a successful detoxification re¬ 
duces the nuiitiier of detoxification op- 
rtunities the next invader must face 
one. In such circumstances the ex¬ 
pected number of invaders that avoid 
random detoxification can be shown to 
be approximately if there 

are Xinvaders and X £ 20, andX-20 
if there are X invaders and X > 20. 

If the probability of a carcinogenic 
respnse is proportional to the number 
of invaders avoid detoxification. 


then the probability of a carcinogenic 
response undergoes an exponmtial in¬ 
crease (nearly linear at very low doses 
and slow overall) as the delivered dose 
of the chemical increases until it gener¬ 
ates approximately 20 invaders. At this 
point further increases in the delivered 
dose correspond to more rapid linear 
increases in the probability of a carcin¬ 
ogenic re^nse. 

In this example the delivered dose 
level correspon^g to tiie generation of 
20 invaders would be called the indi¬ 
viduals susceptibility frontier /. In gen¬ 
eral, if an individual has a nonzero sus¬ 
ceptibility frontier, then/occurs at the 
point where, as the dose increases, the 
effective dose shifts from a slowly in¬ 
creasing function of the delivered dose 
to a more r^idly increasing linear 
function of the delivered dose. 

One form of the Individualized Re¬ 
sponse Model corresponds to a multi- 
nage model, which is linear in the bio¬ 
logically effective dose. 

Mathematically, 

Pit; BEDid.doJ)] “ 1 - exp [-oo- 
Of; X SED(d,4)J)] 

when the susceptibility S is the suscep¬ 
tibility frontier E, and /is the value of 
E. The greek constants oo and a; would 
usually be estimated from the chronic 
dose-response data. 

The functional form of the biologi¬ 
cally effective dose, BED{d,do^, 
would be determined on the buis of 
scientific research ancillary to the re¬ 
sponse frequency tabulation in clinic 
bioassays. Alternatively, the functional 
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How th* tnd<v(dua({z«d 
Rospons* Model works 

In an Individualized Response 
Model the proportion of the 
entire population at risk that is eX' 
pected to develop a specified carcin- 
ogettic response, in the absence of 
competing mortalityi by time t for an 
administered dose d is the weighted 
average of 

Pit; BED(d,do,s)] ■ dose-response 
model in terms of the biologically ef¬ 
fective dose BBD{d,do,s); 

where the weights reflect the propor¬ 
tion of the population that have a 
susceptibility S equal to s and a 
background dose 0 equal to dg. 
Herein. P/t; BED(d,do,s)] denotes 
the probability of the sp«ified re¬ 
sponse occui^ng by time f, in the 
absence of competing mortality for 
the particular biologically effective 
d^ B£D(d,do,s) that arises when 
the administered dose is d, the ran¬ 
dom background dose Dg takes on 
the value dg, and the random sus¬ 
ceptibility S takes on the value s. Fig¬ 
ure 5 provides an overview of the 
components and structure of an Indi¬ 
vidualized Reponse Model. 


biological system through time, P[t; 
BED(d,do,s)] can refiect explicitly not 
only the different stages in the carcino¬ 
genic process, but ^so cell prolifera¬ 
tion, biological processes that are age- 
dependent, and doses that are not 
constant over time. Such models allow 
risk to be characterized in terms that 
reSeci how the probability of a carcino¬ 
genic response changes with time. 

In the multistage theory of carcino¬ 
genesis, a carcinogenic response occurs 
when k successive events stemming 
from a single cell have occurred (d). In 
this context, P[t; BED{d,do,s)] is deter¬ 
mined by the rates at which the number 
of cells completing successive events 
increase over time. Figure 6 provides 
mathematical details for the general 
if-stage carcinogenic process. 

Conceptually, the probability P[t; 
BED(d,do,s)] arises in a very natu^ 
way. For exanqile, if the carcinogenic 
process involves two stages, then the 
process involves normal stem cells, in¬ 
termediate cells that have undergone 
the first of the two carcinogenic events, 
and malignant cells that arise from in¬ 
termediate cells by undergoing the sec¬ 
ond carcinogenic event. Malignant 
cells are assumed to produce the speci¬ 
fied carcinogenic response. 

The number of malignant cells that 
arise at a particular time equals the rate 


at which an intermediate cell gives rise 
to malignant celts at that time multi¬ 
plied by the total number of intermedi¬ 
ate cells alive at that time. The total 
number of intermediate cells at a par¬ 
ticular time, say T, depends primarily 
on the number of intermediate cells that 
arise from normal stem cells at a time 
T’ prior to T and the net growth in that 
number of intermediate cells due to cell 
birth and death between T' and T. 

The number of intermediate cells that 
arises from normal stem cells at a par¬ 
ticular time is the rate at which a nor¬ 
mal stem cell gives rise to intermediate 
cells at that time multiplied by the total 
number of normal stem cells alive at 
that time. 

The entire carcinogenic process and 
its components are all tracked through 
time. That is, the number of cells in 
each stage, the transition rates from 
stage to stage, and the net within-stage 
proliferation rates due to cell birth (rep¬ 
lication) and death or terminal differen¬ 
tiation are evaluated explicitly at each 
point in time. 

Thus the dependence of these factors 
on time and the biologically effective 
dose at that time can be explicitly incor¬ 
porated. Because the biologically effec¬ 
tive dose is calculated at each time, 
time-dependent changes in the adminis¬ 
tered dose, the background dose, and 


form of the biologically effective dose 
is estimated during the fitting of the 
chronic dose-response data. The esti¬ 
mated biologically effective dose could 
increase linearly with the delivered 
dose for delivered doses exceeding the 
susceptibility frontier and increase ex¬ 
ponentially (with associated low-dose 
linearity) for delivered doses less than 
the susceptibility frontier. 

Modeling through time 

The Individualized Response 
ModeTs probability that a randomly se¬ 
lected individual from the population 
will have a specified carcinogenic re¬ 
sponse 1^ time r for an administered 
dose is the expectation of P[t; BED 
id,dg,s)] with respect to the distribution 
of the background dose Dg and the sus¬ 
ceptibility S in the population. The rela¬ 
tionship, P[r; BED(d,do,sy], between 
the likelihood of a response by time t 
and the biologically effective dose 
BED{d,dg,s) for given values of dg for 
Dg and of s for 5, is averaged in the 
Individualized Response Model. This 
model can be a better, more coniplete 
reflection of the available scientific in¬ 
formation on the careinogenic risk as¬ 
sociated with a specific situation if the 
dynamic nature of the biological proc¬ 
esses involved is explicitly reflected in 
P[t; BED(d,do,s)\- 

By modeling the progression of the 


FIGURE 5 

Individualized Response Model: components and strocture 
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the uncertainty in the high* to low-dose 
extrapolation portion of quantitative 
risk assessment. 

By incorporating modeling dirough 
time, the dose-response reli^onship 
for a given biologi^y effective dose 


For example, Thorslund et al. (7) 
suggest a correspondence between 
mo(^ of action in a two-stage (1^=2) 
carcinogenic process and the param¬ 
eters shown in the Figure 6. (See box.) 

Such studies would allow the model 


the susceptibility can also be incorpo¬ 
rated. 

Other models are special cases of the 
form for Fir; BED(d,(io,s)] shown in 
Figure 6 (5-7). These other models as¬ 
sume that the background dose Do and 
the susceptibSi^ 5 have the same val¬ 
ues for every individual and that the 
susceptibili^ frontier F is zero for 
every individual. 

Benefits 

There are several substantial benefits 
to using the Individualized Response 
Model. By explicitly including the 
interindividual variation in susceptibil¬ 
ity and background dose, quantitative 
risk assessment can more accurately re¬ 
flect the population at risk. Particular 
distributions of susceptibility and back¬ 
ground doses within the population can 
be studied and directly incorporated 
into the model. This removes some of 
the unreasonable elements present in 
most quantitative risk assessments, 
namely, the assumption that everyone 
in the population is as sensitive as the 
most sensitive individual. 

The model uses dose on the most rel¬ 
evant scale—the biologically effective 
dose scale—and incorporates changes 
in the dose over time. Physiological 
and pharmacokinetic research can be 
used to provide information on the de¬ 
livered dose corresponding to a particu¬ 
lar route of exposure, the interaction of 
administered and background doses, 
and any effects of susceptibility on the 
delivery process. 

Molecular biology and a careful anal¬ 
ysis of a chemical's mechanistic behav¬ 
ior can establish the nature of intracel¬ 
lular susceptibility and the presence or 
absence of a susceptibili^ frontier. 
Such research can also indicate which 
phenomena are age-dependent and the 
functional form of such dependencies. 
Much of this dosimetry may either be 
already available or obtainable from 
short-term experiments (5). 

Furthermore, the number and range 
of nonzero dow levels for which the 
biologically efiective dose is obtainable 
are generally much greater than the few 
experiment^ dose levels considered in 
a chronic bioassay. Hence, the relation¬ 
ship between the administered dose and 
the biologically efiective dose should 
be obtainable over a wide range of ad- 
mimstered doses, and the biologically 
efiective dose corresponding to the ac¬ 
tual exposure levels of interest can be 
more reasonably determined. 

In addition, tlie research results con¬ 
cerning the likely functional depen¬ 
dence of the dose-response model 
parameters on the dose assume that foe 
dose is expressed on foe biologically 
efiective dose scale. The use of foe bio- 
logicaUy efiective dose should reduce 


can be more biologically motivated. 
Such modeling allows P[t: BED- 
(d,do,s)] for given BED(d,do,s) to be 
analyz^ in biologically meaningfol 
terms, reflecting both foe qualitative 
and quantitative mechanistic informa¬ 
tion generated by foe biologist and foe 
toxicologist. 


structure and Issues such as linearity 
versus nonlinearity to be determined by 
scientific inquiry rather than by as¬ 
sumption or ^licy. The ability to incor¬ 
porate supplementary scientific investi¬ 
gations should promote increased data 
acquisition, provide broader scientific 
pa^cipation in quantitative risk assess- 


10M Env^ran. Se«. ISeAnof.. MV. 2). Mo. >4, 1M7 


Source: https://www.industrydocuments.ucsf.edu/docs/ffkxOOOO 






Corrcspondanca between 
inodes of action* and 
parameters in Figure 6 

• Cocarcinogens, which induce re¬ 
generative hyperplasia as the nor¬ 
mal consequence of a tissue's at¬ 
tempt to repair toxic damage, 
affect Co: 

• Initiators, which induce mutation, 
oncogene activation, or chromo- 
somai translocation, affect Mo', 

• Promoters, which increase the 
number of transformed preneo- 
plastic cells through clonal expan¬ 
sion. affect Gt: 

• Completers, which increase the 
rate of transformation from the pe¬ 
nultimate stage to a malignant 
cell, affect M,-. 

• Inhibitors, which reduce the rate 
of transformation to malignancy 
by different mechanisms, cause 
decreases in Cg, Mg, M,, or Gt- 

Because the model uses biologically 
meaningful parameters, its time- and 
dose-dependence may be investi¬ 
gated in special studies supplemen¬ 
tal to die collection of response fre¬ 
quencies in chronic bioassays. 

•In t iwxug* (*-2) cvcinoganic proctis. 


ment, and instill greater confidence in 
the results. 

Some of the greatest benefits of the 
new modeling occur when the human 
risk characterization involves extn^ 
lalion. Modeling with biologically de¬ 
fined components permits the identified 
differences between species, exposure 
route.s. and exposure time frames 
(chronic, short-term, and only adult 
years) to be explicitly incorporated into 
the model. 

The functional form of the particular 
components where differences do oc¬ 
cur, such as the biologically effective 
dose and the number of normal stem 
cells (Co), can be updated to reflect the 
new situation. Also, the distributions of 
background dose and susceptibility for 
the new situation can replace the old 
distributions. 

Individualized Response Models can 
incorporate new resulu on the individ¬ 
ual components of the dose-response 
model and utilize the specific functional 
forms of these components for particu¬ 
lar situations. Thus these new models 
provide human risk characterizations 
that are more up-to-date, accurate, and 
peirinent to the particular exposure sit¬ 
uations being evaluated. 

Prognosis 

Quantitative cancer dose-response 
modeling is the subject of intensive re¬ 
search. New dose-response models can 
incorporate research on more biologi¬ 


cally relevant dose scales. The new 
models still involve some simplifying 
assumptions and are not exact biologi¬ 
cal models; however, the extensions de¬ 
scribed make it possible to incorporate 
more of the current scientific under¬ 
standing of carcinogenesis and more bi¬ 
ological reality. 

These more complete biologically 
based models provide more accurate 
and informative quantitative risk as¬ 
sessments to both the public and risk 
managers. Furthermore, the new 
models may inspire research directions 
that will yield the experimental data 
needed to more accurately determine 
the dose-response relationship for a 
particular carcinogen. The new capa¬ 
bilities fm* dose-response modeling and 
an expanding data base provide the po¬ 
tential for substantial improvements in 
quantitative cancer risk assessment. 

References 

(1) Harris, C. C. la Human Cardnogenesii; 
Harris, C. C. Auirup, H. N., Eds.; Aca¬ 
demic; New York, 1983: pp. 941-70. 

(2) Cleaver, J. E. In Banbury Rtport 16: Gt- 
ntlic Hiriabillry in Responses to Chemical 
Eiposure; Omenn. G. S.; Gelboin. H. V,, 
Eds.; Cold Spring Harbor: Cold Spring Ha^ 
bor. N.Y., 1984; pp. 225-28. 

(3) Omenn. C. S. In Banbury Report 16: Ge¬ 
netic Ksriabitiry in Responses to Chemical 
Exposure; Omenn. C. S.: Gelboin, H. V., 
Eds.: Cold Spring Harbor: Cold Spring Har¬ 
bor. N.Y.. 1984; pp. 3-13. 

(4) Armiiage. P.; Doll, R. tn Proceedings of 
the Fourth Berkeley Symposium on Mathe¬ 
matical Statistics and Probability; Neyman, 
J., Ed.; University of California wess: 
Berkeley and Los Angeles. 1961; pp. 19-38. 

(3) Moolgavkar, S. H.: Venzon, D. J. Math. 
Biosei. 1979, 47. 55-77. 

(6) Moolgavkar. S. H.; Knudson, A. G. J. 
Nat. Cancer Inst. 1981, 66. 1037-32. 

(7) Thorslund, T. W.; Brown, C. C.; 
Charnleyu. C. J. Risk Anaivsis, 1987, 7, 
109-19. 

(8) Menzel. D. B. Environ. Set. Technol. 
1987,2/, 944-50, 



Robert L Sielken, Jr., is a biosiaiisitcian 
with a Ph.D. in probability and statistics 
from Florida State University In 1985 he 
formed Sielken, Inc, to do fitll-time statisti¬ 
cal research and consulting in the quanti¬ 
tative dose-response modeling and risk 
characterization portions of risk assess¬ 
ment. He serves as a consultant to profes¬ 
sional societies, industry, and state and 
federal governments. Earlier in his career 
Sielken was a teacher and active re¬ 
searcher in the Department of Statistics at 
Texas AAM University. 



TheChemistiyof 
Acid Rain: Sources 
and Atmospheric 
Processes 


T his new book takes a probing look at a 
high-priority environmental problem— 
the sources and chemisny of acidic 
species in the atmosphere. The editors be¬ 
gin this volume with an overview chapter 
that provides a historical perspective and 
sjmmarizes the understanding that has 
been developed over the last decade. You'll 
read subsequent chapters covering all as¬ 
pects of this problem... from the accuracy 
of field measurements to highly sophisti¬ 
cated modeling. Although the problems 
and solutions presented in this book are 
wide-ranging, you’ll find the 27 chapters di¬ 
vided Into seven specific subject areas; 

• CiMnI 
• Rictptw Model] 

• Ckxid Chimistiy ind niyjkt 
• KindiCI 

• Wet tnd Oiy Deposition 
• Experimenai Method] 

• FoiMUmenul eroceuc] 

This volume provides essential information 
for those involved in the acid deposition 
held; atmospheric scientists, ecologists, en¬ 
vironmental scientists, government re¬ 
searchers. regulators, l^lslators. and gen¬ 
eral scientists. 
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